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ABSTRACT 

The baryon acoustic oscillation (BAO) experiment requires a sufficiently dense sam- 
pling of large-scale structure tracers with spectroscopic redshift, which is observa- 
tionally expensive especially at high redshifts z ^ 1. Here we present an alternative 
route of the BAO analysis that uses the cross-correlation of sparse spectroscopic trac- 
ers with a much denser photometric sample, where the spectroscopic tracers can be 
quasars or bright, rare galaxies that are easier to access spectroscopically. We show 
that measurements of the cross-correlation as a function of the transverse comoving 
separation rather than the angular separation avoid a smearing of the BAO feature 
without mixing the different scales at different redshifts in the projection, even for 
a wide redshift slice Az ~ 1. The bias, scatter, and catastrophic redshift errors of 
the photometric sample affect only the overall normalization of the cross-correlation 
which can be marginalized over when constraining the angular diameter distance. As 
a specific example, we forecast an expected accuracy of the BAO geometrical test 
via the cross-correlation of the SDSS and BOSS spectroscopic quasar sample with 
a dense photometric galaxy sample that is assumed to have a full overlap with the 
SDSS/BOSS survey region. We show that this cross-correlation BAO analysis allows 
us to measure the angular diameter distances to a fractional accuracy of about 10% 



at each redshift bin over 1 



3, if the photometric redshift errors of the galaxies, 



fi 2 /(l + z), are better than 10 — 20% level. 
Key words: cosmology: distance scale — large-scale structure of the Universe 



1 INTRODUCTION 



Various cosmological da ta sets such as th e cosmic microwave 
background (CMB; lHinshaw et al. II2012D. the type-la super- 



nova observations (Riessetal 



1998; Schmidt et al. 199 



Perlmutter et all Il999l ; iKessler et all [20091 : ISuzuki et all 



20121 ) a nd the baryon acoust ic oscillation (BAO) exper- 



iments jEisenstein etall 120051; JPercival et all 120071. l20ld ; 



iBlake et al.ll201ll ; lBeutler et alll201ll ; I Anderson et alH2012T ) 

have shown increasing evidence that the cosmic expansion 
today is in the accelerating expansion phase. The cosmic 
acceleration is the most tantalizing problem in cosmology. 

Among others, the BAO experiment is recognized as 
one of the most promising geometrical tests, because it rests 
on the physics of the CMB anisotropics in the early universe, 
which is remarkably well described by the linearized pertur- 
bation theory for a A-dominated cold dark matter model 
(ACDM). The tight coupling between baryons and photons 
prior to the decoupling epoch of z ~ 1100 leaves a character- 
istic imprint on the pattern of large-scale structure tracers 
such as galaxies and quasars - the so-called BAO scale. The 



BAO scale is now preci sely constrained as ~ 150 Mpc from 
the CMB observations |Hinshaw et al.l l2012h which can be 
used as a "standard ruler" to infer the cosmological dis- 
tances from the observed correlation functio n of the tracers 
jSeo fc Eisenstei j|2003l ; iHu fc Haimar]|2003l ) . 

The BAO measurements mostly utilize a large data 
from wide-area redshift surveys of galaxies, such as the 
6dFG^|, the Sloan Digital Sky Survey (SDSSfl the Baryon 
Oscillation Spectroscopic Survey (BOSSfl and the WiggleZ 
survejQ- With the success of these surveys, there are sev- 
eral future BAO surveys targeting higher redshift ranges 
of z >l, including the HETDEX survejfl the Extended 
Baryon Oscillation Spectroscopic Survey (eBOSS the Big- 



http:/ /www. aao.gov.au/6dFGS/ 

http:/ /www. sdss.org/ 

http: / /www. sdss3.org/surveys/boss.php 

http://wigglez.swin.edu.au/site/ 

http: / /hetdex.org 

http: / /www. sdss3.org/future/eboss.php 
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BOSifl the Sub aru Prime Focus Spectrograph projeejfl 
(|Ellis et al.ll2012l) , and the satellite Euclid missions. How- 
ever, extending the BAO experiment to higher redshifts 
(z ~ 2 — 4) is observationally expensive, because the tar- 
get galaxies become increasingly fainter and spectroscopic 
surveys of such faint galaxies covering wide-areas are quite 

time-consuming^- 

In addition to these spectroscopic BAO analysis, there 
have been attempts to measure the BAO feature in 
the correlation function of photometric galaxy samples 
llPadmanabhan et ahll2007l ; iBlake et aT1l2007l ; ICarnero et all 
l2012l ; ISeo et al.ll2012h . A wide-area, multi-colour photomet- 
ric survey is relatively easy to carry out compared to a spec- 
troscopic survey of similar area coverage. In fact, there are 
many planned imaging surveys, including the Subaru Hy- 
per Suprime-Cam (HSC) SurvejO, the Dark Energy Sur- 
vey (DES0, Euclid, and the LSST projecQ for which 
the primary science driver is weak lensing based cosmology. 
However, the photometric BAO experiments are challenging 
for several reasons. First, the photometric BAO analysis is 
based on the angular correlation function of the galaxies, 
which is by nature two-dimensional and loses (wash out) 
the clustering information in the line-of-sight direction that 
also causes a reduction in the correlation function ampli- 
tude. Second, the projection along the line-of-sight smears 
the BAO feature in the angular power spectrum. Third, the 
BAO feature inferred from the photometric samples can be 
significantly affected by statistical and systematic (catas- 
trophic) errors of the photometric redshifts (photo-z's). For 
example, including the photo- z outliers in the analysis can 
easily induce a bias in the BAO peaks, which in turn causes 
a bias in the inferred distance. 

In this paper, we propose to use the cross-correlation 
between the spectroscopic and photometric tracers of large- 
scale structure as an alternative BAO method. This method 
is particularly useful when sampling of the spectroscopic 
tracers is too sparse to measure the BAO feature via its 
auto-correlation analysis. Since a photometric survey usu- 
ally has a much denser sampling, the cross-correlation mit- 
igates the shot noise contamination to improve clustering 
measurements. We argue that smearing due to the line-of- 
sight projection can be avoided by measuring the correlation 
function as a function of the transverse comoving separation 
rather than the angular separation. As a specific example, 
we consider the cross-correlation of the SDSS and BOSS 
spectroscopic sample of quasars with photometric galaxies 
to estimate an expected accuracy of the derivable geomet- 
rical test. The SDSS quasars are bright and can easily be 
observed spectroscopically, but have a too sparse sampling 
for the auto-correlation analysis. When making the forecast, 
we also include the broad-band shape information of the 



7 http://bigboss.lbl.gov/ 

8 http://sumire.ipmu.jp/en/2652 

9 http: / / sci.esa.int / science-e / www/object / index. cfm?fobjectid=48983 

10 We note that the BAO feature at z < 2.3 is recently de- 
tected from the three-dimensional correlation function of the 
Lyman-a forests that are identified i n the BOSS quasar spectra 
llBusca et alj|2012l : ISlosar et al]|2013t) . 

11 http:/ /www. naoj.org/Projccts/HSC/index. html 

12 http:/ /www. darkenergysurvey.org 

13 http://www.lsst.org/lsst/ 



cross-correlation in addition to the BAO feature (also see 
ICoorav et aTll200ll ). 

This paper is organized as follows. In Section [2j we de- 
scribe explicit expressions for the cross-correlation analysis 
as a function of the transverse comoving separation as well as 
its counterpart in Fourier space, and also derive the covari- 
ance matrix. We show our basic results in Section [3] In Sec- 
tion^ we show an expected accuracy of the geometrical test 
via the use of the cross-correlation of the SDSS/BOSS spec- 
troscopic quasar sample with a dense photometric galaxy 
sample. We summarize our results in Section [5] Unless oth- 
erwise stated, we employ a concordance ACDM model, with 
flmoh 2 — 0.137 and Q.^h 2 = 0.023 for the matter and baryon 
physical density parameters, Qa = 0.721 for the cosmologi- 
cal constant assuming a flat geometry, and A s = 2.43 x 10~ 9 , 
n 3 — 0.96, and a 3 = for the primordial power spectrum 
parameters. 



2 BAO FEATURE IN THE PROJECTED 
CORRELATION FUNCTION 

2.1 Transverse cross-correlation function and the 
power spectrum 

In this paper, we consider a method using the cross- 
correlation of a photometric sample with a spectroscopic 
sample for measuring the BAO scale. A key idea is to con- 
sider the cross-correlation measured as a function of the 
transverse comoving separation rather than of an angular 
separation 



w(R) 



1 



[(n s (~f s ;z s )n p (~f p )) - l] 



(1) 



where quantities with subscripts "s" and "p" denote those 
for spectroscopic and photometric samples, respectively; 
n s ( / y s ;z s ) and np(-y p ) are the projected number density 
fields for the spectroscopic and photometric samples in the 
directions of -y s and 7 p on the celestial sphere, respectively; 
z 3 is the redshift of each spectroscopic object. Thus the den- 
sity field of the spectroscopic sample is described as a func- 
tion of both z s and the angular position. The transverse 
radius R is defined in terms of their observed angular posi- 
tions of ~y s and ~f p and the redshift z 3 as 



7? = d A (z s ) cos (■y a ■ 7 p ) ~ d A (z s )\6s 



(2) 



The quantity <1a{zs) is the comoving angular diameter dis- 
tance to each spectroscopic object. Note that this conver- 
sion requires to assume a background cosmological model. 
The unit vector on the celestial sphere, 7, is given as 
7 = (sini? cos <p, sin sin cp, cos ■&). In the last equality on 
the r.h.s. of Eq. we used the flat-sky approximatiorF 4 !. 
Observationally, the cross-correlation is estimated from the 
average of all the pairs separated by the same separation R 
within a given width, compared to the cross-correlation of 
the spectroscopic sample with random catalogues that are 



14 In the flat-sky approximation, the unit vector is expanded 
around the North pole as 7 ~ (1? cos <p, 1? sin ip, 1), and the 
two-dimensional flat-space vector can be defined as 6 = 
(1? cos (p, 1? sin ip). 



© 0000 RAS, MNRAS 000, 000-000 



BAO with the cross-correlation 3 



constructed based on the same selection function as in the 
photometric catalogue. 

A notable advantage of the i?-average over the angle 
average is that it can preserve the physical scales inherent 
in large-scale structure such as the scale of baryon acous- 
tic oscillation in which we are interested, if the assumed 
cosmological model for the R-6 conversion is closer to the 
true underlying cosmology. On the other hand, the ^-average 
mixes different scales in large-scale structure, thus smearing 
the BAO scale in the observed cross-correlation function. 
We emphasize that this _R-average is useful when a spectro- 
scopic catalogue is available for the cross-correlation mea- 
surement. In contrast, in the case of the auto-correlation 
analysis of photometric samples, the conversion from 9 to R 
is severely affected by photo- z uncertainties, which can lead 
to a smearing or systematic offset of the BAO feature. 

We can express the projected cross-correlation function 
in terms of the power spectrum as follows. First, in the case 
that the spectroscopic sample has a redshift distribution, 
the projected cross-correlation can be expressed as 



simple form: 



)(R) = / dz s ps(z s ) w(6; z 3 )\ R=d , )e . 
Jo 



(3) 



where w(8;z 3 ) is the angular cross-correlation function be- 
tween a spectroscopic sample at redshift z s , Ps(z s ) is the 
redshift distribution of the spectroscopic sample, normal- 
ized as J °° dz 3 Ps(z 3 ) — 1, and the average with the no- 



tation \ R - dA ( z \0 indicates that the redshift average for a 
given R is done by averaging the angular correlation func- 
tion w(9; z s ) under the condition R = (Ia(zs)O according to 
the discussion around Eqs. (QJ and p]l. w(9; z s ) is defined in 
terms of the angular power spectrum as 



w(8; z s 



-t-Y(2l + l)C sp {l;z s )P,{cos6). (4) 



Here Pi (x) is the i-th order Legendre polynomials. For a flat 
universe, the angular power spectrum is given in terms of the 
three -dimensional po wer spectrum, in a standard manner 
(e.g.. lDodels"onll2003l ). as 

C 3 p(l;z s ) = ~ J drW p (r) J k 2 dk P sp (k; z s , z)ji(kr s )ji(kr), 

(5) 

where P Bp (k;z s ,z) is the three-dimensional cross-power 
spectrum between the spectroscopic objects at redshift z a 
and photometric objects at z; r is the radial distance given 
as a function of redshift for a given cosmology, r = r(z), and 
r s — r(z s ); ji(x) is the l-ih order spherical Bessel function; 
W p (r) is the selection function of the photometric sample, 
normalized as JdrW p (r) = 1 (see below for an example). 
Note that r(z) = dA(z) for a flat universe. 

Using th e flat-sky ap proximation and the Limber's ap- 
proximation (|Limberlll954T ). the projected cross-correlation 
function can be simplified as 

C krlk 

w(R) = I ^Lc sp {k)J {kR), (6) 

where Jo(x) is the zeroth-order Bessel function, and the 
transverse comoving separation separation R is defined for 
spectroscopic redshift z 3 of each sample used in the average. 
The projected cross-power spectrum C ap (k) is given by a 



dz . 



C ap (k) = J dr P s(z) — W p {r)P sp (k;z). (7) 

Note that the power spectrum C ap (k) has a dimension of 
Mpc 2 so that k 2 C sp (k) becomes dimension- less. We use 
C sp (l) for the usual angular power spectrum and C sp {k) for 
the Fourier counterpart of the w(R) throughout the paper. 

We assume that we can, based on photo-z technique, 
select photometric objects that have similar photo-z to the 
spectroscopic redshift. Even in the presence of large photo- 
z errors, the cross-correlation method is very powerful in 
the sense that it can statistically select photometric objects 
that are physically clustering with the sp ectroscopic sample 
(|Newmanll2008l ; iMcQuinn fc Whitdl2013t ). Including photo- 
z bias and outliers in the sample simply dilutes the cross- 
correlation signals, but does not change the shape so that the 
BAO scale is not shifted. If the spectroscopic sample used 
in the cross-correlation measurement is in a narrow range of 
redshifts, [z s ,z 3 + Az s ], the projected power spectrum reads 



C sp {k) 



drW p {r) 



Ar. 



■Psp(k;z s ), (8) 



where r s and r s + Ar s are the radial distances to red- 
shifts z a and z s + Az s , respectively, and we have used 
Ar s p s (zs){dz/dr) = 1. The factor j rs+Ars drW p (r) is a 
fraction of photometric objects among the whole photo- 
metric sample that reside in the spectroscopic redshift bin 
[z s , z s + Az s ] and thus are physically correlated with the 
spectroscopic sample. Therefore the factor gives a dilution 
factor of the cross-correlation signal due to the photo-z er- 
rors. The factor 1/Axs in front of P sp (k) accounts for the 
fact that the cross-correlation amplitude is reduced with in- 
creasing the width of spectroscopic redshift bin. The above 
equation explicitly shows that an inclusion of the photo-z 
outliers does not change the shape of the cross-correlation, 
but simply affects the overall normalization. Also impor- 
tantly, the projected cross-correlation, measured against R 
instead of the angular separation, can measure the three- 
dimensional power spectrum P sp {k; z s ) at spectroscopic red- 
shift and at a particular k. Put another way, the projected 
cross-correlation does not mix the power spectrum of dif- 
ferent Fourier modes, which is not the case for the angu- 
lar power spectrum. It should also be noted that the pro- 
jected cross-correlation is not affected by redshift-space dis- 
tortion due to the peculiar motions of the tracers. In par- 
ticular, the nonlinear redshift-space distortion, the so-called 
Finger-of-God effec t, is very difficult to accurately model 
(Hika ge"et al,ll2012T ). and therefore the cross-correlation may 
have a practical advantage. 

Next let us consider a case that the spectroscopic sam- 
ple is in a given width of the redshift bin, z a = [z s OVJ , z" p ], 
where z a ovl and z" p are the lower and upper bound of the 
redshift bin respectively. Here for simplicity we consider a 
uniform distribution of the spectroscopic sample within the 
given redshift bin 



(z) = 4, Az s 1 







otherwise, 



(9) 



where Az s = — z s ow . The uniform redshift distribution 
is not a critical assumption, and can be easily generalized to 
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Figure 1. The thin solid line shows a toy model of a photometric 
galaxy distribution. Among them, we select galaxies which fall in 
the spectroscopic redshift (z s ) bin. The thick solid lines show the 
underlying true redshift distribution of the photometric samples 
which are defined by the photo-z bins of z p £ [0.6, 1.0], [1.0, 1.8], 
and [1.8,3.2], assuming a photo-z accuracy of cr z = 0.01(1 + z s ) 
(see Eq. HOP . The dotted lines are similar distributions but the 
photo-z errors are degraded to <r z = 0.1(1 + z a ). 



a case that the spectroscopic sample has a redshift distribu- 
tion. For photometric objects used for the cross-correlation 
measurement, we would select the objects if the best-fit 
photo-z's are in the range of the spectroscopic redshift bin. 
Here we employ the simplified assumption that the proba- 
bility for photo-z's obeys a Gaussian distribution 



p p (z\z p ) = 



2ira : 



■ exp 



2a 2 . 



(10) 



where we assumed the photometric objects with the best-fit 
photo-z, z p , obeys a single population, ct z is the lcr photo- 
z error, and z is its true redshift. The probability satisfies 
the normalization dz p p p (z\z p ) = dz p p (z\z p ) — 1. 
Given the distribution, if the photometric objects whose 
photo-z's are in the range of the spectroscopic redshift range, 
Zp G [zs° w , Zg P ], the probability distribution for the true red- 
shift is computed as 



p P (z\z p € [z'° w ,z" p ] 



dz p pp(z\z p ) 



~ [erf(* up ) 



r / 1< 

ert(x 



(ii) 



where erf (a;) is the error function and a; up,/low = (zs P//low — 
z)/y/2az. Taking into account the overall redshift distribu- 
tion of the photometric sample, we can derive the redshift 
distribution of the photometric sample based on the photo-z 
selection: 



3 (z) = in p (z)[erf( 2; up )~erf( a; 10 ")], (12) 

where n p (z) is the redshift distribution of the photometric 
sample for which we a ssume n v (z) = (z 2 /2z q) exp[— z/zp] 
parametrized by zo (see lOguri fc Ta kada 201 1]). Throughout 
the paper we assume zo = 0.4, yielding the mean redshift 
(z) = 1.2, and assume the bias parameter of the photometric 
sample to b p = 1.5. Hence, the selection function of the pho- 
tometric sample in each z s bin used for the cross-correlation 




Figure 2. The histogram shows the redshift distribution of 
the s pectroscopic quasar sample of SPSS DR7 ( Schneide r et alj 
l20ldl and BOSS DR9 jParis et al J 120121) . The BOSS sample 
is multiplied by 3 as the DR9 has completed only the 1/3 of 
the target area. Dashed line shows the assumed photometric 
galaxy distribution, n p (z) defined below Eq. J 1 21) . where we as- 
sumed n p = 50 arcmin" 2 for the total mean number density and 
(z) = 1.2 for the mean redshift. 



is given as 



W p (r) = 



n pSZs (z) dz 



dr ' 



(13) 



where n p< z Zg is the normalization factor, defined as 
n p g Zs = J °°dz n p6Zs (z), so as to satisfy the condition 
/ °° dr W p (r) = 1. We will study how an accuracy of the 
BAO measurement changes with the quantities such as <j z 
and Az s as well as the number densities of spectroscopic and 
photometric samples. 

In Fig. [TJ we show an example of the photometric galaxy 
distributions with different photo-z uncertainties, when the 
photometric galaxies are divided in a redshift bins, which is 
chosen to match with the spectroscopic sample. In this ex- 
ample, we divide the whole sample into three subsamples as 
Zp e [0.6,1.0], [1.0,1.8], [1.8,3.2]. The photo-z errors cause 
a leakage of the photometric galaxies from the spectroscopic 
redshift bin. 



2.2 Covariance matrix 

The error covariance matrix quantifies an accuracy of mea- 
suring the projected cross-correlation for a given survey. 
Since angular scales at different redshifts are scaled to match 
the transverse comoving scale for an assumed cosmologi- 
cal model, the measured projected cross-correlation is two- 
dimensional, given as a function of the comoving scales in 
units of Mpc. Assuming a Gaussian error for the projected 
power spectrum, which is a good approximation at BAO 
scales l|Takahashi et al.1l2009l ). we can extend the standard 
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formula for the covariance matrix for angular power spectra 
|Knoxll995l ) to obtain the covariance matrix of the projected 
cross-correlation function as 



Cov[C 3p (k),C ap (k')] = 



°kk> 



N mode (k) 



C ap {k) 2 + I Css(k) + — ) \C„{k) + 



,(14) 



where JV mo de(fc) is the number of independent Fourier mode 
discriminated by the given survey area defined as 

27rfcAfc 



A, 



»(*) 



2tt 



cU(4 ow )e s 

= 2fcAfc rfA(z'° W ) 2 /sky, 



(15) 



with / s ky being the sky coverage defined as / s k y = f2 s /47r. 
The quantities n 3 and n pSZa are the projected number den- 
sities of the spectroscopic sample and the photometric sam- 
ple having photo-z's values within the spectroscopic red- 
shift bin, respectively. The number densities are in units of 
Mpc -2 . In the above equations, we assumed that the funda- 
mental model of the two-dimensional Fourier decomposition 
is defined as the projected scale at the lowest redshift for a 
given survey area, kf = 27t/[g(a(zs OW )0s]. 

The expression for the covariance matrix (Eq. I14p 
can be used to understand why the cross-correlation can 
be useful for the BAO analysis when the spectroscopic 
catalogue has too sparse sampling of the targets, i.e., 
Css >C l/n s - We assume that the photometric sample has 
a high number density in the spectroscopic redshift bin, 
i.e., C pp 3> l/n p g Za , at the BAO scale, and the cross- 
correlation coefficient between the spectroscopic and pho- 
tometric samples, r = C sp / \JC S sC pp , is of order unity. In 
this case, the covariance for the cross-correlation becomes 
Cov[C sp ,C sp ] oc Cp P /n 3 , and thus the signal-to-noise ra- 
tio is (S/N)2 P = CjpCov -1 oc r 2 C aa n a . This should be 
compared with the case of the auto-power spectra of the 



spectroscopic sample, (S/N) M oc (C a . 



yielding the 



S/N ratio (S/N)'/(S/N)i 



r 2 /(C aa n a ) > 1. Thus the 



cross-correlation can give a higher signal-to-noise ratio for 
a measurement of the projected power spectrum at the 
BAO scales. In practice, however, the spectroscopic sam- 
ple allows a measurement of the three-dimensional power 
spectrum, which contains more Fourier modes than in the 
projected power spectrum. In the next section, we present 
a more quantitative comparison between the methods us- 
ing the cross-correlation and the three-dimensional auto- 
correlation. 

Fig.[2]shows the photometric and spectroscopic samples 
for which we think the cross-correlation method discussed in 
this paper is useful, if the two survey regions are overlapped. 
The spectroscopic quasar catalogues of the SDSS and BOSS 
surveys have a wide coverage of redshift up to z a ~ 4, but 
have a much lower number density than in the photometric 
galaxies available from the upcoming imaging surveys such 
as the Subaru HSC Survey or Euclid. The redshift distri- 
bution of the photometric sample shown in Fig. [2] is deeper 
than what is usually assumed for the HSC survey or Euclid, 
but we note that for the cross-correlation analysis we can 
use fainter galaxies than galaxies used for the weak lensing 
analysis. We also note that our results are not very sensi- 



tive to the choice of the number density distribution of the 
photometric sample. 



3 RESULT 

3.1 Projected power spectra 

In Fig. Owe compare auto- and cross-power spectra for spec- 
troscopic and photometric samples at mean redshift z = 1.4. 
Here the cross-power spectra is computed as a function of 
the transverse comoving separation as described in the pre- 
vious section. Here, we consider the redshift bin around 
z = 1.4, z=[z- Az/2, z + Az/2] with widths Az = 0.365 
and 1.05, which correspond to the radial distance widths of 
Ar = 0.5 and 1.5 Mpc//i, respectively. We assumed that we 
can select objects in the redshift bin, based on the spectro- 
scopic redshift or photo-z's. To model the photo- z errors, we 
use the parametrization given in lMa et al.l {2006) as 



A z (l + z s ), 



(16) 



and consider the two cases of A z = 0.05 and 0.3. The BAO 
feature is smeared for the angular auto-power spectra of pho- 
tometric samples, while the BAO feature persists in the pro- 
jected auto- or cross-power spectra using the spectroscopic 
sample (see also below). Comparing the solid and dashed 
curves shows that the larger photo-z errors cause a more sig- 
nificant dilution of the power spectrum amplitudes, thereby 
smearing the BAO oscillatory feature. The figure also shows 
the shot noise levels. For ongoing or upcoming imaging sur- 
veys, we typically have more than 10 4 deg -2 galaxies (see 
Fig. [2), and thus the power spectrum measurement has a 
sufficient signal to noise ratio. However, the photo-z errors 
dilute the spectrum amplitude and smear the BAO feature, 
suggesting that it would be difficult to use the angular auto- 
power spectrum of the photometric galaxies for an unbiased 
BAO geometrical test, as we will discuss below. While the 
projected cross-power spectra shown in the figure also have 
a more diluted amplitude with increasing photo-z uncertain- 
ties (see Eq. [BJ, we can still use the unsmeared BAO feature 
for estimating the angular diameter distance. 



3.2 Prospect for the cross-correlation BAO 
measurement 

In this section, we study prospects for the use of the pro- 
jected power spectrum for measuring the BAO feature. In 
Fig. 21 we show the projected cross-power spectrum as 
a function of the transverse wavenumber, divided by the 
no-wiggle power spectrum (with the BAO feature being 
smoothed out), in order to highligh t the BAO feature. Note 
that we used the transfer function in lEisenstein fc Hul (|l998h 
to compute the no-wiggle spectrum for the same cosmologi- 
cal model. Although we assume a linear bias multiplicative- 
factor for both the spectroscopic and photometric samples, 
we include the effect of nonlinear clustering on the matter 
power spectrum, us ing the publicly available code, RegPT 
(|Taruva et al.l l2012h . that includes up to the 2-loop order 
contributions based on the refined perturbation theory. We 
show the cross-power spectra up to a certain maximum 
wavenumber, fc max , which is determined so that the nonlin- 
ear matter power spectrum at the mean redshift is expected 
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10 100 1000 10 100 1000 10 100 1000 



multipole 1 



Figure 3. Comparison of the projected and angular power spectra at the mean redshift z = 1.4. {Left): The angular auto-power spectrum 
of photometric samples in the photo-z bins z v = [1.4 — Az/2, 1.4 + Az/2], The thick and thin curves show the power spectra for the bin 
width of Az p ~ 0.365 and 1.05, which corresponds to the width of the comoving radial distance, Ar = 0.5 and 1.5 Gpc/h, respectively. 
The solid and dashed curves are the spectra assuming the photo-z accuracies of <r z /(l + z) = 0.05 or 0.3, receptively. Each thin dotted 
lines show the shot noise level for the photometric samples, which typically have the projected number density more than 10 4 deg~ 2 for 
an imaging survey we are interested in. {Middle): Similar to the left panel, but for the cross-power spectra between the spectroscopic 
and photometric samples, as a function of the transverse comoving separation (Eq. [7J, where the transverse mode k is rescaled to the 
multipole via the distance to the spectroscopic sample by I = kr{z = 1.4) for an illustrative purpose. The solid and dashed curves are 
for the photo-z accuracies of the photometric galaxies, as in the left panel. The cross-correlation preserves the BAO wiggles compared to 
the left panel. {Right): The projected auto-power spectra for the spectroscopic samples. The figure shows that, for a spectroscopic survey 
with a small number density n a < 10 2 deg - 2 , the BAO wiggles in the auto-spectrum are difficult to measure due to the significant shot 
noise. 



to be accurate to within a 1%-le vel accuracy in the ampl i- 
tude compared to the simulation l|Taruva et al . 2009, 2012). 
The figure clearly shows that the projected cross-power spec- 
trum preserves the BAO feature, even for a wide redshift bin 
of the projection. On the other hand, the BAO feature is 
smeared in the angular correlation. We also notice that, for 
a higher redshift slice, the BAO feature remains up to the 
greater wavenumber due to the less evolving nonlinearities. 

We estimate prospects for detecting the BAO feature in 
the projected cross-spectrum by using the ^-difference be- 
tween the power spectra with and without the BAO feature 



[C S p{ki) — C"^(fci)] 
Cov [C sp (ki), Csp(ki)] ' 



(17) 



where C sp and C™ are the cross-power spectra with and 
without the BAO feature, and the summation is up to the 
maximum wavenumber determined as in Fig. [4] Note that 
A\ 2 does not include the broad-band shape information 
of the cross-power spectrum, and quantifies a significance 
of detecting the BAO feature in the cross-power spectrum, 
assuming that the spectrum with the BAO wiggles is the 
underlying true spectrum (see the next section for a more 
quantitative forecast of the BAO analysis). The denomina- 
tor is the covariance matrix (Eq. I14p for which we assumed 
the Gaussian error. To compute A\ 2 , we assume that a 
spectroscopic sample has the projected number density of 
n s — (20Az) deg _2 in the redshift bin, where Az denotes the 
redshift bin width. For the photometric sample, we employ 
n* ot = 1.8 x 10 5 deg -2 (50 arcmin -2 ) for the total number 



density. These numbers roughly resemble the SDSS/BOSS 
quasar spectroscopic sample and the Subaru HSC Survey 
or Euclid imaging surveys, respectively. We here assume a 
full-sky coverage (/sky = 1) for both the spectroscopic and 
photometric catalogues; A\ 2 oc ,f s ky We then assume that 
we can select the photometric objects, which have photo-z's 
in the spectroscopic redshift bin, and take into account the 
redshift distribution of photometric objects as well as the 
effect of photo- z errors using the method in Sec. 12.11 



Fig.[S]shows the A\ 2 values for the cross-power spectra 
assuming a various combination of the survey parameters. 
If the two surveys have a sufficiently wide area coverage 
for their overlapping region, the projected cross-power spec- 
trum allows a detection of the BAO feature. We compare 
the results with a BAO analysis for the spectroscopic sam- 
ple alone. Similarly to Eq. (|17|l , we can define the differential 
X 2 to quantify the sensitivity of the three-dimensional power 
spectrum to the BAO feature 



2nk 2 dk 
(2tt) 3 



p 8 {k)-pr(k) 



P s {k) + n~ 3 



(18) 



where P s {k) is the three-dimensional power spectrum of the 
spectroscopic sample and Kurvoy is the survey volu me. We 
ignor e the redshift-space distortion for simplicity l|Kaiserl 
Il987l ). The figure shows that, if the photo-z accuracies of 
a{z p )/(l + z) are better than 10 — 20%, the cross-correlation 
can achieve a more significant detection of the BAO feature 
than in the three-dimensional power spectrum. 
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Figure 4. The projected powe r spectrum divided by the no-wiggle linear power spectrum in order to highlight the BAO feature, where 
we used lEisenstein fc Hul jl998) to compute the no-wiggle spectrum. We consider z„ = 2.5 (top panel), 1.4 (middle), and 0.8 (bottom), 
respectively, for the mean redshift of the projection. The thick curves are the spectra computed taking account of the non-linearity of 
the matter power spectrum (see text for the details), while thin curves show the cross-power spectra in linear theory. Note that we 
assume the linear bias parameter of b = 1 for both the spectroscopic and photometric samples for simplicity. The spectra are plotted up 
to fc max , where the nonlinear power spectra are expected to be accurate to within 1% at each mean redshift compared to simulations. 
For comparison, we also show the angular cross-power spectra projected over the redshift slice of Ar = 0.9 Gpc/h (short- dashed) and 
Ar = 1.5 Gpc/h (long-dashed) around each mean redshift, respectively, which are plotted against the wavenumber using the conversion 
kr(z) = I. The dotted curves are the nonlinear power spectrum using the no-wiggle linear power spectrum for the input spectrum. 



4 GEOMETRICAL TEST WITH THE 
CROSS-CORRELATION METHOD 

In this section, we present more quantitative estimates on 
the power of the cross-correlation method for determin- 
ing the angular diameter distance. For this forecast, in 



contrast to the preceding section, we include the broad- 
band shape informatio n of the cross-power spe ctrum, ex- 
tending the method in ISeo fc Eisensteinl (|2003h to a two- 
dimensional cross-correlation analysis. As a specific exam- 
ple, here we consider the cross-correlation BAO analysis 
assuming the SDSS/BOSS spectroscopic quasar catalogues 
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Figure 5. An expected significance of the BAO detection, A\ 2 fEq, H7jl . for the cross-correlation analysis with different combinations of 
spectroscopic and photometric samples. We estimate the significance by comparing the cross-power spectra with and without the BAO 
wiggles, as in Fig. [3] but do not include the broad-band shape of the power spectra. In each panel, the thick solid curves show the Ax 2 
values for the projected cross-power spectrum (C sp (k)), the thick dashed curves are for the angular cross-power spectrum (C sp (l)), and 
the thin horizontal line is for the projected auto-correlation of the spectroscopic sample (C ss (k)). For comparison, we also show the result 
when the BAO feature is extracted from the three-dimensional power spectrum analysis (P:jj)(k)), which is estimated using Eq. ]18t . 
The number density of the spectroscopic sample is fixed to (20 Az) deg -2 , and the total number density of the photometric sample is 
assumed to n* ot = 50 arcmin - 2 . Results are shown for three mean redshift, z g = 2.5 (top panels), 1.4 (middle), and 0.8 (bottom). Left: 
A\ 2 is calculated under the conditions that the photo- z accuracy is fixed to A z =0.1 and the redshift bin width is varied from Ar = 0.1 
to 1.5 Gpc//i. Right: A\ 2 is calculated with a fixed redshift bin width Ar = 1.5 Gpc/h, but with varying photo-z accuracies. 



l|Schneider et al.ll2010l ; iParis et al.ll2012r ) for a spectroscopic 
sample, as shown in Fig. For a photometric sample, we 
assume a hypothetical imaging sample which has a full over- 
lap with the spectroscopic survey region; we assume the to- 
tal area of 10,000 deg 2 for the overlapping area. We consider 
six redshift bins with the mean redshifts ranging from 0.7 
to 2.9. The projected number density in each bin is esti- 
mated using the redshift distribution in Fig. [5] We use the 
bias parameters of the quasars in each re dshift bin based 
on the measurement by iRoss et al.l l|2009l ). For the photo- 
metric sample, we again assume the total number density 
of n* ot = 50 arcmin -2 , and compute the number density in 
each redshift bin taking into account the photo-z error (see 
Sec. 12. ip . Table [1] summarizes the set of the survey param- 
eters. 

The cross-correlation is measured as a function of the 



transverse separation between the pairs of the spectroscopic 
and photometric objects. The transverse separation, the sep- 
aration distance between each pair perpendicular to the line- 
of-sight direction, can be inferred from the observed angular 
separation on the sky, R oc A9 (see Eq. [5J. For this con- 
version, we need to assume a reference cosmological model 
to relate the observable A8 to the quantity R. Thus the 
transverse wavenumber is given as 

k ±M =-^Lk±. (19) 

The quantities with "ref" are the quantities from the ob- 
servable assuming a "reference" cosmological model, and 
the quantities without the subscript denote the underlying 
true quantities. Since the reference cosmological model as- 
sumed generally differs from the underlying true cosmology, 
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Table 1. A summary of survey parameters we consider for the forecast, and the expected fractional errors of determining the angular 
diameter distance, cr(Dj\)/ Da, including marginalization over the other parameters. Here we consider the SDSS/BOSS spectroscopic 
quasar catalogue for the spectroscopic sample, and the Subaru HSC- or Euclid-type galaxy sample for the photometric sample. z s and Az a 
are the mean redshift and the redshift width for each redshift bin of the spectroscopic sample taken in the hypothetical cross-correlation 
analysis. 6 S , /3 and n a are the linear bias, the linear redshift-space distortion and the number density in each redshift bin (see text for the 
details). fc max is the maximum wavelength used for the Fisher matrix analysis. For each redshift bin, we cross-correlate the spectroscopic 
sample with the photometric galaxies based on their photo-z's assuming the photo-2 errors of A z = 0.01, 0.1 and 0.3, respectively (see 
Eg. 1161 ■ n p is the number density of the photometric galaxies in each redshift bin (see Eg. 1121 1. The last column (&b a /Da) shows the 
expected error on the angular diameter distance measurement in each bin. For comparison, we also show the expected error when using 
the three-dimensional auto-power spectrum of the spectroscopic sample ("spec auto-correlation"). 



it causes an apparent distortion in the c ross-power spectrum 
- the so-called Alcock-Paczynski test (lAlcock fe Paczvnskil 
1 19791 ; lEisenstein fe Hulll99sl : ISeo fe Eisensteinll2003l ). Thus 
the observed cross-power spectrum is given as 

C? p s {k^ ci - z) = D t! xe f*t C sp {k; z) + V s {z), (20) 

where V s {z) is the residual shot noisj^l. Since we consider 
a wide redshift bin for the projection of the spectroscopic 
sample, we can safely ignore the redshift-space distortion 
effect. 

To make a parameter forecast, we include the following 
set of parameters: 

6 = {Q. m , n m /i 2 , Q. b h 2 , A s ,n a , a a , D A (zi)> A{zi), V s (zi)}, 

(21) 

15 Suppose that the spectroscopic and photometric samples re- 
side in their host haloes, which have the number densities of n^i 
and n^2, and assume that some fractions of the two samples have 
the common host haloes which have the number density of n c . In 
this case, the shot noise for the cross-power spectrum is found to 
be proportional to n c j (n^in^)- 



where £l m h 2 and Q.bh 2 are the matter and baryon density 
parameters today, A a , is the amplitude of the primordial 
curvature perturbation at fc p i v = 0.005 Mpc, and n 3 and a a 
are the tilt and runn ing of the primordial power spectrum 
ijKomatsu et al.ll201ll ). The parameter Da{z%) is the angular 
diameter distance to the i-th redshift bin which is treated as 
an independen t parameter from other cosmological par am- 
eters (see e.g.. ISeo fe Eisensteinl [iool ; lEllis et al.ll20ll for 
details). We also include the normalization parameter A(zi) 
which models an uncertainty in the normalization of the 
cross-power spectrum in each redshift bin due to unknown 
bias uncertainties for both the photometric (b p ) and spec- 
troscopic (b s ) samples. As discussed above, the marginaliza- 
tion over A(zi) also takes account of photo-z uncertainties. 
Vs(zi) is a nuisance parameter to model the residual shot 
noise parameter. In addition to these parameters, we include 
the optical depth and angular diameter distance to the last 
scattering surface, r and Da,cmb, respectively to the Fisher 
matrix to describe the CMB prior. 

The full Fisher matrix can be expressed by a simple 
sum of two Fisher matrices, F = F CMB + F GC , where F cc 
denotes the Fisher matrix from the cross-correlation mea- 



© 0000 RAS, MNRAS 000, 000-000 



10 A. J. Nishizawa et al. 



surement: 



f: 



a/3 



E E 



9^; 



Cov" 1 



W B 



(22) 



where Cov is the covariance matrix given by Eq. (I14p . The 
maximum wavenumber fc max is set to the maximum scale up 
to which the nonlinear matter power spectrum at the mean 
redshift is expected to be accurate to within 1% level as in 
Fig.U 

For comparison, we also show a forecast for us- 
ing the three-dimensional power spectrum of the spectro- 
scopic sample to estimate the cosmologica l distances, H(zj) 
and D A(zj). W e follow the met hods in ISeo fc Eisensteinl 
|2007ft (also see lEllis et alj|2012h. We model the redshift- 



Kaised Il987t) and its nonlinear 
20071 ) for the three-dimensional 
additional parameters; = 



space d istortion (RSD; 
effects l|Eisenstein et al] 
power spectrum with 
d In D(zi)/d In a/b 3 and H(zi). The fiducial value of /3 is 
listed in Table[T] However the results are not sensitive to the 
details, because the power spectrum information at relevant 
wavenumber bins is limited by the shot noise contamination 
for the sparse spectroscopic sample we are interested in. 

Table [T] and Fig. [6] show an expected accuracy of the 
angular diameter distance measurement in each redshift 
bin via the cross-correlation method. The cross-correlation 
method allows an improvement in the geometrical test com- 
pared to the three-dimensional auto-power spectrum anal- 
ysis, by reducing the shot noise contamination. For the 
SDSS/BOSS spectroscopic quasar catalogues, the cross- 
correlation method enables the fractional accuracy better 
than 10% in each redshift bin, if the photometric galaxy 
survey has a full overlap with the SDSS/BOSS footprints 
and we can select adequate galaxy samples whose photo-z 
errors are better than A z = 0.1 — 0.2 (see Eg. I16|) . Also, an 
advantage of this method is to determine the angular diame- 
ter distance up to a high redshift of z ~ 3, where the cosmic 
expansion is well in the decelerating expansion phase. 



5 SUMMARY 

In this paper, we have studied how the cross-correlation be- 
tween the spectroscopic and photometric samples can be 
used for the two-dimensional BAO experiment. We have 
shown that, with the aid of the spectroscopic redshift, the 
cross-correlation preserves the BAO feature in the probed 
transverse scales, even for the projection over different red- 
shifts such as Az ~ 1, while the angular (cross)-correlation 
suffers from a smearing of the BAO feature due to unavoid- 
able photo-z errors that cause a mixing of the different phys- 
ical scales in a particular angular scale (see Fig. [3]). There 
are several notable advantages of this method. First, the 
cross-correlation allows a significant reduce in the shot noise 
contamination to the measurement. Second, any statistical 
or systematic (catastrophic) photo-z errors affect only the 
overall normalization of the cross-correlation function, and 
do not change the shape of the power spectrum. 

The cross-correlation method can be useful, if the spec- 
troscopic sample has a wide coverage of redshift, but does 
not have a sufficiently high number density for the BAO 
experiment via the auto-correlation analysis. As a specific 
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Figure 6. Expected accuracies on the angular diameter dis- 
tance measurements with the cross-correlation BAO analysis as 
in Table[T] We show the results for the photo-z accuracies of both 
A z = 0.01 and 0.1, denoted as C sp (fc) (0.01) or C sp (fc)(0.1), respec- 
tively. The outermost boxes show the expected accuracies when 
the auto-correlation power spectrum of the spectroscopic quasar 
catalogue is used. The spectroscopic sample is divided into 6 sub- 
samples with their spectroscopic redshifts. The solid and dashed 
curves show the changes in the angular diameter distance when 
the dark energy equation of state (w) is changed to —2 and —1/3 
from w = — 1 (the cosmological constant). 



example, we have considered the SDSS and BOSS spectro- 
scopic quasar sample to estimate the feasibility of the cross- 
correlation method, motivated by the fact that wide-area 
imaging surveys such as the Subaru HSC Survey and Eu- 
clid are being planned, which overlap with the SDSS survey 
footprints (see Fig. [2}. Here the SDSS quasar sample has a 
wide redshift coverage of < z 4 and wide area coverage 
of about 10000 deg 2 , but has too small number density of 
~ 10 2 deg -2 per unit redshift interval to implement the BAO 
experiment via the auto-correlation analysis. On the other 
hand, the planned imaging surveys likely provide a much 
denser sampling of galaxies such as 10 3 — 10 5 deg -2 over the 
redshift range. We have shown that the cross-correlation al- 
lows a more accurate BAO measurement over 0.7 < z £ 3 
than in the auto-correlation of the spectroscopic sample or 
the angular power spectrum of the photometric galaxies (see 
Figs. [5] and [6] and Table[T}, if the photometric redshift is rea- 
sonably good, 10 — 20% level in the fractional accuracy, in 
order not to have a severe dilution in the measured cross- 
correlation. As shown in Fig. [6] the better photo-z accuracy 
of a z /(l + z s ) — 1% does not improve constraints on Da sig- 
nificantly compared to the 10% photo-z accuracy. Hence the 
10 — 20% of the photo-z accuracy is sufficient for the cross- 
correlation BAO study, which can easily be achieved for 
the current and upcoming multi-band photometric galaxy 
surveys. Large-scale structure is more in the linear regime 
at higher redshift due to less evolving nonlinearities, and 
therefore we can use the shape of the cross-power spec- 
trum, including the BAO feature, up to the higher maximum 
wavenumber, thereby yielding a more constraining power of 
the angular diameter distance. 

The expected accuracy of the angular distance mea- 
surement in Fig. [6] is from both the BAO feature and the 
broad shape of the power spectrum. The projected cross- 
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correlation allows us to measure the shape of the three- 
dimensional power spectrum, although the overall normal- 
ization is affected by photo- z errors. Hence, the method can 
also be used to constrain the tile and running index of the 
primordial power spectrum. Also, as an ultimate possibil- 
ity, the cross-correlation method may enable to use the ob- 
served radius of dark matter haloes in the projected dis- 
tance. If we have a good knowledge on the virial radius of 
dark matter haloes as well as have a good estimator of halo 
masses, observed the virial radius can be used to infer the 
angular diameter distance. This is relevant for cluster-shear 
weak lensin g, which probes the ha lo and dark matter cross- 
correlation (|Oguri fc Takadall2011r ). Given that the clusters 
have follow-up spectroscopic redshifts, we can expect a high- 
precision measurement of the halo-matter cross correlation 
at small scales down to a few Mpc, which correspond to 
the virial radii of massive haloes. Thus the virial radius may 
serve as another standard ruler that can lead to even higher- 
precision measurements of the angular diameter distances. 
This is an interesting possibility and may worth exploring 
further. 

We note that the cross-correlation technique developed 
in this paper can also be used to constrain the primor- 
dial non-Ganssianity via measurement s of the largest-scale 
cross-power spectra jDalal et al. 20081 ; iTaruva et al.l 120081 ; 
iMcDonaldl |200S| ; ISlosar et al.ll200Sft . Again, by measuring 
the cross-correlation as a function of the transverse comov- 
ing separation, we can avoid the the smearing effect due 
to the projection which reduces the enhanced power at the 
largest scales of k £j 0.01/i/Mpc. 

Spectroscopic observations of quasars, or more gener- 
ally bright, rare galaxies, are relatively inexpensive in terms 
of the observation time needed for a given telescope. Such 
objects are also very interesting subjects for astronomical 
studies. The method developed in this paper can add a cos- 
mological science case when combined with wide-area imag- 
ing surveys that have an overlap with the spectroscopic sur- 
vey. The method is useful when designing joint spectroscopic 
and photometric surveys including a science case of the two- 
dimensional BAO analysis. 
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